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*---  '  to  Levy  [7]  and  Sanathanan  and  Koerner  [8]  has  been  suc- 


A  Dominant-Data  Matching  Method  for  Digital  Control 
Systems  Modeling  and  Design 

L.  S.  SHIEH,  Y.  F.  CHANG,  AND  R.  E.  YATES 

Abstract— A  d— l Iwt-dsti  matching  method  Is  presented  for 
obtaining  a  red  nerd -order  paiae- transfer  Auction  from  either  a  high- 
order  contlnnnne  data  transfer  Auction  or  a  high-order  paiae— 
transfer  Auction,  and  far  Identifying  the  poise-transfer  reaction  of  a 
system  from  available  experimental  time  and  fttqmcy  response 
data.  The  method  may  also  he  applied  to  the  digital  control  systems 
design  problem  with  various  sampling  periods. 


I.  INTRODUCTION 

Most  practical  industrial  circuits  and  control  systems  are 
continuous-time  systems  for  which  analog  filters  and  control¬ 
lers  are  employed  to  improve  performance.  The  recent  avail¬ 
ability  of  high-performance  low-cost  microprocessors  and  as¬ 
sociated  electronics  has  led  to  replacement  of  many  continu¬ 
ous  systems  with  system  employing  digital  filters  and  con¬ 
trollers.  Many  techniques  have  been  developed  for  digital  con¬ 
trol  systems  design  (ll-[4).  Among  them  the  w-domain  bi¬ 
linear  transformation  is  often  applied  to  design  industrial 
digital  controllers.  However,  this  method  is  graphical  and  in¬ 
volves  “cut-and-try”  procedures.  Recently,  Kuo  1 5  J  and 
others  developed  an  optimal  discrete-time  data  matching 
method  for  the  redesign  of  a  continuous-data  system.  Con¬ 
stant  controllers  instead  of  dynamic  digital  controllers  are 
mainly  employed  in  these  designs.  As  a  result,  good  perform¬ 
ances  of  redesigned  systems  can  be  achieved  if  the  frequency 
of  the  input  signal  is  sufficiently  lower  than  the  sampling  fre¬ 
quency.  As  an  alternate  to  Kuo’s  time-domain  approach, 
Rattan  and  Yeh  [6]  have  given  an  elegant  frequency-domain 
method  for  the  redesign  of  continuous-data  systems.  The 
method  of  weighted  least  squares  complex-curve  fitting  due 
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cessfully  extended  in  the  2-domain  to  determine  a  dynamic 
digital  controller.  As  a  result  of  these  efforts,  better  perform¬ 
ance  of  redesigned  systems  can  be  achieved.  In  this  corre¬ 
spondence,  a  computer-aided  method  is  proposed  for  match¬ 
ing  the  dominant  data  of  a  high-order  continuous-data  sys¬ 
tem,  or  a  discrete-data  system,  with  the  dominant  response 
of  a  low-order  digital  system  replacement.  Also,  methods  are 
given  for  system  identification  and  digital  controller  design  of 
these  systems. 

II.  DOMINANT-DATA  MATCHING  METHOD 

The  characteristics  of  a  control  system  or  a  filter  are  often 
expressed  by  either  a  time-response  curve,  or  a  frequency-re¬ 
sponse  curve  or  a  set  of  poles  and  zeros  in  the  complex  plane, 
or  both.  The  quantitative  description  of  transient  behavior  is 
represented  by  its  time-domain  control  specification  [9]  (for 
example,  the  percentage  overshoot  and  the  rise  time)  and  by 
the  frequency-domain  control  specifications  (for  example,  the 
maximum  value  of  the  closed-loop  frequency  response  and  the 
bandwidth).  These  specifications  which  are  defined  for  analog 
control  systems  [9]  can  be  considered  as  specifications  of 
digital  control  systems.  This  is  because  a  digital  control  system 
can  be  viewed  as  a  continuous-data  system  in  the  frequency 
domain  when  2  =  e'^T  where  T  is  a  sampling  period. 

Some  empirical  observations  or  rules  of  thumb  due  to 
Axelby  [10]  and  Truxal  [11]  that  link  the  specifications  of 
the  continuous-data  systems  in  both  the  time  and  frequency 
domains  can  be  found  in  [  10] ,  [  1 1 J . 

Recently  Shieh  et  al.  [12]  have  studied  the  relationships 
between  the  complex-domain  specifications  (the  damping 
ratio  and  the  undamped  natural  angular  frequency)  and  the 
time-domain  and  frequency-domain  specifications.  Based  on 
the  results  f  12] ,  the  time  and  complex-domain  control  speci¬ 
fications  can  be  converted  into  the  equivalent  frequency- 
domain  specifications  and  vice  versa.  These  control  specifica¬ 
tions  are  used  as  dominant  data  of  a  continuous  or  discrete- 
data  control  system. 

Our  new  dominant-data  matching  method  matches  the 
dominant  data  of  a  continuous-data  or  a  discrete-data  system 
to  those  of  the  newly  designed  or  modeled  discrete-data  sys¬ 
tem.  The  steps  involved  are  as  follows: 

Step  1:  Determine  a  set  of  dominant  frequency-response 
data  from  the  assigned  or  obtained  time-domain  and  complex- 
domain  specifications  by  using  the  rules  and  results  in  [10]- 
[12]. 

Step  2:  Assume  a  fixed  configuration  digital  system  and 
controllers  with  unknown  constants.  Determine  the  open-loop 
and  the  overall  pulse-transfer  function  of  the  system. 

Step  3:  Formulate  a  set  of  linear/nonlinear  equations  by 
matching  the  unknown  constants  of  the  pulse-transfer  func¬ 
tion  and  the  assigned  dominant  data.  Solve  the  equations  by 
using  the  multidimensional  Newton-Raphson  method  [13], 
available  as  a  library  computer  program  package  (called  the 
Z  system)  in  many  digital  computers  [  14] . 

Step  4:  Estimate  initial  value  for  the  numerical  solution  of 
the  Newton-Raphson  method  by  constructing  a  crude  pulse- 
transfer  function  which  can  be  obtained  by  a  complex-curve 
fitting  method. 

Step  5:  Compare  the  results  with  the  assigned  specifications. 
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Pig.  1.  Block  diagram  of  a  missile  pitch  control  system. 


*k  »  M  *1  *1  fcfc  *| 


"S|  |fe  k|  ||  fek  -  |Y  |fe  |k  |5t 

Fig.  2.  The  Nyquist  plots  of  various  example  open-loop  transfer  functions. 


III.  MODELING  A  REDUCED-ORDER  PULSE-TRANSFER 
FUNCTION 

We  use  a  real  stabilized  pitch  control  system  of  a  semiac¬ 
tive  terminal  homing  missile  [IS]  as  an  illustrative  model  to 
show  that  the  characteristics  of  the  transient-state  response  of 
a  system  can  be  estimated  from  the  dominant  frequency-re¬ 
sponse  data  and  the  applications  of  the  proposed  method  to 
the  identification  and  model  reduction  problems.  A  block 
diagram  of  the  missile  system  is  shown  in  Fig.  1 .  The  closed- 
loop  high-order  transfer  function  is 


G0(s)  =  transfer  function  of  the  actuator  and  aerodynamics 
of  the  missile  system 

-  324  332.316(s  +  0.1933X*  +  6SX*  +  1S00) 

"  s(s-  2.92 1 X* + 3  •  1 75X* + 87.9  ±/95 ,5X*  +  H  2.5X*  + 1 385  ) 

(lc) 


Hg(s)  =  transfer  function  of  the  gyro  =  1 


T(s)  Gc(sX?q(s)  Ge{s) 

R(s)  l+Gc(syG oltyijs)  1  +  Ge(s)  ~ 

where 


(la) 


(Id) 

G,(s)  =  Gc(s)Go(j)  =  unstable  open-loop  transfer  function  of  the 
existing  stabilized  system. 


Gc(s)  =  stabilization  filter 
1.6 


(g+'XiH 

Wf 


(le) 

The  Nyquist  plots  of  Ge(s)  and  G0(r)  are  shown  in  Fi*.  2.  The 
dominant  data  of  G,(a)  are: 

1 )  Real  and  imaginary  parts  of  G,(s)  a 1 1  =  /a)  =  /0  are 


Ob) 


Re|G,yO)]  =*-2.103  817 
Im  (0,00)1  =  ». 


(2«) 
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2)  Gain  margin  G,m  of  this  system  Ge(ju„)  is  type  “I”  systems  and  we  need  to  match  the  dominant  data 

of  Ge(ju )  as  to  =  0  in  (2a)  and  the  dominant  data  of  G/z) 
^  ]  1  (  |  1  1  |  as  to  =  0  or  z  =  e'uT  ~  1  in  (3b)  or  (4),  the  expression  for 

Gem  =  r  a, ,  \  =  b-  ir  a, ,  m  ~  17*7  Gr(z)  in  <3b) is  modified  as  follows: 

!G,0^)l  Re  [G,0«,)]  1.5|  Substituting  z  =  z*  +  l  into  (3b)  yields 


where  the  phase-cross-over  frequency  to,  is  given  by 


C0jr  —  1.9  rad/s,  such  that  Z-G,(/cj,  )  =  - 1 80° . 

(2c) 

The  equivalent  real  and  imaginary  parts  of  G,(/co,)  at 
1 .9  rad/s  are 

We  = 

where 

Re  [G^O'w,,)]  =-1.507  944 

(2d) 

lm  [G,0«„)]  =  -0.006  490  205  . 

Xn * 

(2e)  p 

3)  Phase  margin  <t>em  of  the  system  G,(/w)  is 

<t>em  =  180°  +  LGe(juc)  a  5.7787° 

(20 

x.  *  +  x„  ,  *z*  +  -  +  x 0*z*p 

Cr(z*)  =  — £ - — - - - — 

z*(yq*  +yq- 

=  e_iZ*~l  +  e0  +ejz*  +  e2z*J  +  — 


P 

=  2  *<• 

i=0 


*P  i* =2  “p-i'  y<>* = 2^'- 
*=  i  <~o 


3V.*  =  2'V. 

»=  1 

,  ,  ,  ,  f-1  =Xp*lyq*  ande0  =(.vg*x„_  i *xp*)lyq**,fHc. 

where  the  gain-cross-over  frequency  toc  is  given  by  toc  s  3.2 

rad/s  so  that  Equating  the  respective  real  and  imaginary  parts  of  G,(z) 

for  to  =  0  and  those  of  G,(z*)  for  to  =  0  gives 

|G,0wc)l  =  l-  (2g) 

Re  [Gr(z)l  |I= ,  +/0  =  Re  (Gr(z*)J  ,  =/0  =  e0  (5b) 

Equivalent  real  and  imaginary  parts  of  G,(/6Jc )  are 

and 

Re  [G,0'toc)l  =  -0.993  914  3  (2h) 

lm  (Gr(z)l  l*=  i  +/0  =lm  [Gr(z*)]  |,.=I_l=/0  =  (5c) 

1m  [G,</«c)l  =  -0.095  474  78.  (2i) 

Equations  (5b)  and  (5c)  imply  that  e0  in  (5b)  is  the  asymp- 
It  is  required  to  determine  a  reduced-order  pulse-transfer  func-  totic  line  of  the  type  “1”  systems  at  low  frequencies, 
tion  such  that  the  characteristics  of  the  identified  discrete-data  In  the  frequency  domain,  (4)  can  be  expressed  in  an  alter- 
model  agree  as  closely  as  possible  with  those  of  the  high-order  native  form  as  follows, 
continuous-data  system.  Let  us  define 

Let  the  required  overall  pulse-transfer  function  be 


7V(z)  = - 

1  +Gr(z) 

Grfuk,vk)  =  — — 


z  =  e/w*r  =  cos  Co* T  +  j  sin  uikT~  uk  +  jvk 
and  substituting  z  =  uk  +  jvk  into  (4),  we  have 


(xquS-xqvS  +  x,u„  +  x2)+/(2x0ukvk  +xtvk) 


t("k  —  lX«k2  — +.V|«*  +y2)-vk{2uicvk+yiuk)\+muk-lX2ukvk+ylvk)+vk(uk1  -vk2  +y,uk+y2)) 


=  Rk+  Vk 


where  the  open-loop  pulse-transfer  function  Grfz)  is 


Gr(z)  = 


x0zp+xizp  1  +  - +xp-,z +  Xp 
(Z- lX-Po*’  +  ■+J’e-i«+3'e) 


G^z)  is  assigned  to  be  a  type  “1”  system  because  G,(z)  in 
(le)  is  a  type  "1”  system.  To  match  the  five  dominant  data 
in  (2)  we  choose  q  =  2,p  =  2,  andy0  =  1.  Thus  G,(z)  becomes 

XnZ1  +  XiZ  +X, 

SX*)-  (4) 


(*-  IX*  +y  jz  +y2) 


where  x;  and  yi  in  (4)  are  unknown  constants  to  be  deter¬ 
mined.  The  goal  is  to  determine  the  unknown  constants  xt  and 
yt  in  Grfz)  so  that  G,(z)  as  z  =  e^wT  matches  the  dominant 
data  in  (2).  The  sampling  period  T  (=0.008  s)  and  w,  (=250tr 
rad/s)  are  chosen  to  be  synchronized  with  the  1 25-Hz  pulse- 
width  modulated  actuator  [16].  Since  Ge(s)  and  G,fz)  are 


where  cuk  are  specific  frequencies  and  Rk  4  Re  [G^u*,!!*)] , 
/ic  =  lm  [Gr(uk,  u*)| .  If  Rk  and  Ik  are  the  known  or  assigned 
values  at  frequencies  «*,  we  can  obtain  two  linear  equations. 
First,  we  multiply  both  sides  of  (6b)  by  the  common  de¬ 
nominator,  then  we  separate  the  real  and  imaginary  parts  and 
then  equate  the  respective  real  and  imaginary  parts.  Thus  we 
have 

/K*o.*i.*2.3'i.2’z) 

=  (*oUk*  -*o P*1  +x,u*  +*j) 

-  iXn*J  -  t>k2  +y2) 

~vk(2ukvk  +y,u*))  +/*[(«*  -  lX2u*u* 

+  y,u»)  +  M«**  +yj)l 

=  0  (6ci 
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ounuT 


■V. 


» 


and 


-Or 0«*o*  +x,vK) 

-**!(«*  -  1X2u*d*  +>i«*)  +  d*(»<*3  ~vk3 
+  y,»k  +>-})l  -/*!(««*  ~  IX«*J  -w*1 

+  J-IR*  *yi)-ok(2Mkuk 

-0.  (6d) 

Using  the  expreask  in  (3b),  (6c),  and  (6d)  and  the  assigned 
dominant  data  in  \  (2d),  (2e).  (2h),  and  (2i),  we  can 
formulate  one  nonlii  equation  and  four  linear  equations 
ff.xlty,)~  0  for*  “  1. 

The  Newton-Raphscn  i..ethod  available  as  a  library  compu¬ 
ter  program  in  moat  digital  computers  ( 14)  can  be  applied  to 
solve  these  nonlinear  equations.  However,  as  is  well  known, 
the  Newton-Raphson  method  will  converge  to  a  desired  solu¬ 
tion  for  a  small  range  of  starting  values  or  initial  solution  esti¬ 
mates.  To  improve  the  convergence  and  to  obtain  the  set  of 
desired  solutions,  we  offer  the  following  method  for  initial 
estimates. 

Since  y()  *  0  is  nonlinear  and  ft(xi,  yj)  -  0,  i  =  2, 
— ,  5  are  linear  equations,  we  linearize  f\(xi,y()  -  0  by  choos¬ 
ing  a  very  low  frequency.  For  example,  if  we  choose  u*  = 
0.01  6  ug  oii  then  Rk  &  Ro.oi  *  —2.1,  /ft  A  /q.oi  =  40.1 73 
19,  si*  A  u0  0,  =  cos  <00.017*  =  0.999  999  50,  and  vk  A 
Do. oi  *  sin  Wo.oi T  *  8  X  I0~*.  Solving  /j (or/,  yi)  ■  0  and 
fAxt,  Yt)  m  0,  /  *  2,  ,  5  for  the  unknown  constants  x{  (de¬ 

fined  as  x/*)  and  yt  (defined  as  y/*)  we  have  the  initial  esti¬ 
mates  x/*  and  yi*.  Using  these  values  as  initial  estimates  for 
the  solutions  of  /X*i.  >i)  *  0  and  using  the  Newton-Raphson 
method  we  obtain  the  exact  solution  xt  and  yt  at  the  second 
iteration  with  error  tolerance  of  10~*.  The  desired  open-loop 
pulse-transfer  function  is 


0.006  792  596z*  -  0.012  335  992*  +  0.005  545  31 1  4 

C,(|)a - ; - r-  . .  — . .  . . 

*3  -  2.998  541  2* 3  +  2.9964864*  -  0.007945  2 

(7) 


A  Nyquist  plot  of  G^z)  is  shown  in  Fig.  2.  The  plot  matches 
closely  that  of  Ge(r)  not  only  at  the  dominant  frequencies, 
but  also  at  others.  The  GAz)  is  seen  to  be  a  good  reduced- 
order  digital  model  of  the  original  1  lth-order  unstable  analog 
system  Ge(s).  This  is  the  contribution  of  our  new  method  be¬ 
cause  there  are  no  known  effective  continuous-discrete  model 
conversion  and  model  reduction  methods  for  unstable  sys¬ 
tems.  The  resulting  closed-loop  pulse-transfer  function  which 
is  the  reduced-order  discrete-data  model  of  the  original  high- 
order  continuous  data  system  is 


7*)  = 


GAz) 

1  +GAz) 


0.006  792  596 z3  -  0.012  335  992z  +  0.005  545  31 1  4 
z3  -  2.991  748  524Z3  +  2.984  150  408z  -  0.992  399  888  6 


(8) 


Since  the  assigned  dominant  data  are  the  steady-state  fre¬ 
quency  response,  it  is  interesting  to  compare  responses  of  the 
1  lth-order  continuous  function  Tc(r)  in  (la)  and  the  3rd- 
order  discrete  function  TAz)  in  (8)  shown  in  Fig.  3.  Observe 
that  both  the  transient  response  and  steady-state  response  of 
the  reduced-order  model  TAz)  are  excellent  matches  of  the 
original  high-order  system.  This  indicates  that  the  dynamic 
characteristics  of  the  system  (for  example,  peak  value  time 
and  overshoot,  which  may  not  occur  at  the  sampling  time)  are 
indirectly  controlled  by  the  assignment  of  the  gain-cross-over 
frequency  and  the  phase  margin.  This  is  a  major  advantage  of 
our  new  method.  Also  note  that  the  reduced-order  model 
gives  an  excellent  approximation  of  the  original  system  when 
driven  by  high-frequency  input  signals. 

To  determine  the  initial  estimates  xf  and  yt*,  a  general 
formulation  of  a  set  of  linear  equations  can  be  constructed 
from  the  following  complex-curve  fitting  method. 

Consider  the  pulse-transfer  function 


G{t) 


x0*zm  +xl»zm-‘  +■••  +  *„' 

_  -j. 


V^Tz*  +  V, 


,n~i 


+  -  +  R . 


(9a) 


V 
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where  >0*  =  1  and  xt*  and  >/*  are  unknown  constants  to  be 
determined.  Substituting  zr  -  eiroi*T  =  cos  rukT+j  sin  r<jk  T 
into  (9a)  gives 

G(e"w*r) 

«  m 

2  xi*  cos  (m  —  /)«* T  +  ;‘2  *i*  sin  (m  —  t)ukT 

_  fco _ i=o _ 

—  »  n 

2)  y,*  cos  (n-[)u>kT  +  /'  2  yf  sin  ( n-t)u}kT 

1=0  1=0 

=  R(co*)  +  //(co*)  =  R*  +  //*  (9b) 

where  R *  and  /*  are  the  real  and  imaginary  parts  of  the  trans¬ 
fer  function  at  the  experimental  frequencies  or  the  assigned 
frequencies  oo*.  After  multiplying  both  sides  of  (9b)  by  the 
common  denominator  and  separating  the  real  and  imaginary 
parts,  we  equate  the  respective  real  and  imaginary  parts.  This 
yields  the  following  matrix  equation: 


Fig.  4.  Block  diagram  of  a  digital  pitch  control  system. 


where 

G*C0(r)  =  (l-r-')Z^iGo(s)j 

Gel(z)£Gc(z)CkG6(z) 

I  —e~,T 

Gk(s)  = -  =  the  zeroth-order  hold 

s 

T  =  the  sampling  period  =  0.008  s. 


cos moo, T cos (m  —  l)oo, T  ••  1  (— Rt  cos(n  —  l)to,T  +  /|  sin(n  —  l)o>i7^  ••  (—, R,  cosoo,r  +  /,  sin  w,r) 
sin  moo,  r sin  (m  —  l)co,  T  ••  0  (— R,  sin  (n  —  l)o>i  T— /,  cos  (n  —  l)co,  T)  ••  ( — /i a  sin  co,  T  — /,  cos  to,  T) 


cos mo),T cos (m  —  1  Jcj,/1  ••  1  (-Rfcos(«  -  l^T  +  Z/Sin  (n  —  l)to,r)  •  •  (~Rt  cos  oj,T  + 1, sin  oo,r) 
sin  mu/T sin  (m  —  l)u/T  •*  0  (-R,  sin(n  -  l)cj,7'  — /, cos(n  —  1  )to,r)  ••  (-R, sin  cos  co/7’) 


(/to  cos  nco0r— 10  sin  ncj0T) 

e 

e 

e 

(R0  sin  nwo T  + /0  cos  noi^T) 

-K< 

v  • 

• 

• 

~h 

yi* 

(Rt  cos  nco/7- — /,  sin  mo,T) 

• 

yz* 

(Rt  sin  nutT  +  It  cos  nu,T) 

• 

• 

- 

(9c) 


Substituting  the  selected  (n  +  m  +  1 )  frequency-response  data 
into  (9c),  we  can  solve  for  the  required  (n+m  +  1)  unknown 
constants  xt*  and  yt*. 

IV.  DIGITAL  CONTROL  SYSTEM  DESIGN 


Gc(z)  is  the  desired  digital  controller  and  is 


Gc(z)  = 


x0z2  +x,z+x2 
z2  +  y,z+y2 


(11) 


Consider  the  pitch  control  transfer  function  of  the  missile 
system  of  (la).  The  unity-feedback  system  without  the  stabili¬ 
zation  filter  Gc(j)  is  unstable,  and  a  rate  gyro  is  not  available 
for  this  example  system.  It  is  required  to  design  a  digital  con¬ 
troller  Gc(z)  instead  of  an  analog  controller  Gc(i)  such  that 
the  designed  system  has  the  exact  control  speciflcations  [9] 
of  the  original  stabilized  continuous-data  system  given  in  (2). 
Furthermore,  the  response  G«(/u)  at  co  -  140  rad/s  =  co140  is 
chosen  as  a  dominant  data  constraint  because  the  system  has 
an  inherent  high  frequency  signal  component  at  co14g,  This  is 
a  digital  redesign  problem.  The  structure  of  the  digital  control 
system  is  shown  in  Fig.  4.  The  closed-loop  pulse-transfer  func¬ 
tion  of  the  desired  digital  system  becomes 


Y(r)  Ge(z)GkG0(z)  y  C«,(z) 

R(z)M  I  +  Gc(z)GkGo(z)  +  G,t(i) 


(10) 


where  x;  and  y,  are  unknown  constants  to  be  determined.  Be¬ 
cause  Gc(z)  is  a  forward  controller,  the  equations  ffat%  yt)  - 
0  can  be  formulated  from  the  following: 


Gc(*)l«=»/w*r 


_ j gut 

CsG0(r)ll.r/wtr 


(12) 


where  co0  =  0,  w,.9  =  1.9,  *  3.2,  and  cot40  *=  140. 

Notice  that  G4G0(r)  for  s  =  e*“*r  is  not  equal  to  G0(a) 
for  s  =  /co*  unless  T  =  0.  Using  the  dominant  data  of  (2a), 
(2h),  and  (2i),  the  required  response  at  U|40,  or  Re 
lGc(e^w,40~)J  =  26.951  878  and  fan  {Ge(e'w»4or))  * 
19.196  865,  and  the  relationships  expreseed  in  (5b)  and  (12) 
yield  a  set  of  equations  ffat,  yj)  *  0  for  1  *  1, 2,  — ,  5. 

The  set  of  linear  and  nonlinear  equations  can  be  solved 
using  the  Newton-Raphson  method.  The  initial  estimates 
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for  the  Newton-Raphson  solution  may  be  determined  from 
(9c).  Another  linear  equation  fi*(xi,  yd  =  0,  instead  of 
/iC*;,  yd  —  0  can  be  constructed  to  yield  five  linear  equa¬ 
tions  with  five  unknown  constants  (x/*  and  yj*).  Ge(jw)  at 
co  —  0.01  4  «o  o  j  is  used  in  this  case.  Substituting  Re 
ICc(e/"oo»7')|  =  1.596  1120  and  Im  [Gc^O-OiT1)]  = 
6.409  642  into  (12)  we  get  the  linear  equation.  Using  xt* 
and  yi*  obtained  from  (9c)  and  the  Newton-Raphson 
method,  we  obtain  the  solution  x /  and  yt  of  /K*/,  yd  =  0  at 
the  fifth  iteration  with  the  error  tolerance  of  10-6.  The  re¬ 
quired  digital  compensator  is 


newly  designed  digital  controller  Gc*(z)  is 
Gc*(z) 

13.170  704z3  -25.531  43Qz2  +  14  629  635z  -  2.268  545  1 
z3  -0.374  248  41z2  -0.327  570  47z  -  0.297  947  56 

(16) 

The  closed-loop  pulse-transfer  function  is 


7V2(*)  = 


G,2(z) 

1  +C,2(z) 


(17) 


Gf(z)  = 


1 1.869  083z2  -  13.492  37z  + 3.058  4008 
z2  -  0.750  552  99z  +  0.646  992  37 


_ 1 1 .869  083(z  -  0.824  080  67Xz  -  0.3 1 2  685  33) 

(z  -  0.37S  276  50  +/0.71 1  449  17Xz  -  0.375  276  50  — /0.71 1  449  17) 


03) 


The  Nyquist  plot  of  Get(z)  4  Gc(z)G„G0(z)  is  shown  in  Fig. 
2.  It  closely  matches  the  Nyquist  plot  of  Ge(s).  The  unit-step 
responses  of  the  existing  stabilized  continuous-data  system 
Te(s )  in  (la)  and  Te,(z)  in  (10)  are  shown  for  comparison  in 
Fig.  3.  The  time-response  of  the  newly  designed  sampled-data 
system  is  very  close  to  the  existing  stabilized  system.  It  is  in¬ 
teresting  to  note  that  Gc(i )  of  (lb)  is  a  fourth-order  analog 
controller  whereas  the  Gc(z)  of  (13)  is  a  second-order  digital 
controller. 

In  a  large  control  system  it  is  often  difficult  to  select  a  min¬ 
imum  common  sampling  period  among  the  various  subsystems. 
For  example,  the  missile  inner  loop  stability  system  with 
sampling  period  T  ~  0.008  s  is  used  with  a  terminal  guidance 
system.  The  terminal  guidance  system  is  low  pass.  Therefore, 
a  larger  sampling  period  may  be  economically  used  in  this 
system.  If  we  assign  a  larger  sampling  period  Tt  (=0.015  s) 
for  the  outer  guidance  loop,  and  we  desire  a  single  sample 
period,  we  must  raise  the  sampling  period  T  (=0.008  s)  of  the 
actuator  and  inner  loop  from  T  (=0.008)  to  T\  (=0.015). 
Notice  that  the  new  sampling  frequency  w,i  (=2n/Tl  = 
418.88)  >  2u>|4o  (=280  rad/s).  The  modified  open-loop 
pulse  transfer  function  with  T t  =  0.015  s  is 


G„G0*(z)4Z 


(14) 


Since  a  larger  sampling  period  T\  is  used,  we  select  a  third- 
order  digital  controller  Gc*(z)  rather  than  the  second-order 
digital  controller.  Gc(z)  of  (13)  is 


where 


G,2(z)  =  Gt*(z)G*G0*(z). 


The  Nyquist  plot  of  Ge2(z)  shown  in  Fig.  2  matches  very  well 
that  of  Ge(s).  The  unit-step  response  for  Y(z)  in  (17)  is  shown 
in  Fig.  3.  The  time  response  of  reJ(z)  very  closely  matches 
that  of  the  original  system  Te(s).  The  resulting  design  is  seen 
to  be  quite  satisfactory. 

V.  CONCLUSION 

A  dominant-data  matching  method  has  been  given  for  fitting 
the  coefficients  of  a  pulse-transfer  function  from  available 
time  and  frequency  response  data  or  from  assigned  design 
goals  expressed  by  a  set  of  control  specifications.  When  the 
dominant  data  are  obtained  from  a  high-order  continuous-data 
as  well  as  a  discrete-data  system,  our  new  method  has  been 
used  to  determine  a  reduced-order  discrete-data  system.  If  the 
data  are  experimental  time  and  frequency  response  data  of  a 
system  to  be  identified,  our  method  may  be  used  to  identify 
the  pulse-transfer  function.  Also,  the  method  has  been  used 
for  redesigns  of  a  continuous  system  using  a  digital  filter  with 
various  sampling  periods.  The  pulse-transfer  function  ob¬ 
tained  by  our  new  method  has  the  exact  dominant  perform¬ 
ance  of  the  original  or  desired  system.  We  feel  that  the  flexi¬ 
bility  and  accuracy  of  our  new  method  will  have  significant 
practical  advantages  for  the  design  of  digital  control  systems 
and  for  the  implementation  of  a  low-order  digital  controller 
by  using  low  cost  microprocessors. 


Gf*(z) 


x0z 3  +X |Z*  +  x2z  +x3 
z3  +  y,z2  +y2z+y  3 


(IS) 


The  x/  and  y2  are  seven  unknown  constants  to  be  determined. 
G,(/cj) at  u>  =  0  4  w  =  1.9  =  Wj.9,  w  =  3.2  4  cjj.j,  and 
u  =  140  4  W|40  shown  in  (2)  are  used  as  the  dominant  data 
to  determine  xt  and  yt.  Using  the  above  design  procedure,  we 
can  determine  a  set  of  equations  yi)  =  0,  ( =  1 ,  2,  7. 

These  equations  can  be  solved  by  using  the  Newton-Raphson 
method,  with  the  set  of  initial  estimates  obtained  from  (9c), 
The  data  obtained  from  (12)  at  u  =  0.01  4  Uq.oi,  wi.p, 
wj.a,  and  Wjao  ere  used  in  (9c)  to  determine  the  initial  esti¬ 
mates  x/*  and  yi*.  Using  these  values  as  initial  values  for  the 
Newton-Raphson  method,  gives  the  desired  constants  xj  and 
yt  at  the  17th  iteration  with  error  tolerance  of  10~*.  The 
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